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■ Abstract 

■ We show that the model of quantum computation based on density matrices and superop- 
erators can be decomposed in a pure classical (functional) part and an effectful part modeling 
probabilities and measurement. The effectful part can be modeled using a generalization of 
monads called arrows. We express the resulting executable model of quantum computing in the 
programming language Haskell using its special syntax for arrow computations. The embed- 
ding in Haskell is however not perfect: a faithful model of quantum computing requires type 
capabilities which are not directly expressible in Haskell. 



o 

1 Introduction 



A newcomer to the field of quantum computing is immediately overwhelmed with many apparent 
differences with classical computing that suggest that quantum computing might require radically 
new semantic models and programming languages. In some sense this is true for two reasons: 
(1) quantum computing is based on a kind of parallelism caused by the non-local wave character 
of quantum information which is qualitatively different from the classical notion of parallelism, 
and (2) quantum computing has a peculiar notion of observation in which the observed part of 
the quantum state and every other part that is entangled with it immediately lose their wave 
character. Interestingly none of the other differences that are often cited between quantum and 
classical computing are actually relevant semantically. For example, even though we do not often 
think of classical computation as "reversible," it is just as reversible as quantum computing. Both 
can be compiled or explained in terms of reversible circuits jz] , but in neither model should the user 
be required to reason about reversibility. 

The two properties of quantum computing discussed above certainly go beyond "pure" clas- 
sical programming and it has been suspected earlier that they might correspond to some notion 
of computational effect. Following Moggi's influential paper [5], computational effects like assign- 
ments, exceptions, non-determinism, etc. could all be modeled using the categorical construction of 
a monad. This construction has been internalized in the programming language Haskell as a tool to 
elegantly express computational effects within the context of a pure functional language. Since the 
work of Moggi, several natural notions of computational effects were discovered which could only 
be expressed as generalizations of monads. Of particular importance to us, is the generalization of 
monads known as arrows [7j which is also internalized in the programming language Haskell. 

In an early paper, Mu and Bird (2001) showed that quantum parallelism is almost a monad. We 
expand and build on this observation as follows. First the traditional model of quantum computing 
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cannot even express measurements, so we use a known more general model using density matrices 
and superoperators. After expressing this model in Haskell, we establish that the superoperators 
used to express all quantum computations and measurements are indeed an instance of the concept 
of arrows (with a small caveat). In particular the construction clarifies the crucial need for some 
form of linear typing: arrow computations must be required to use every quantum value or else the 
computations produce results that arc inconsistent with quantum mechanics. 

In summary, our construction relates "exotic" quantum features to well-understood semantic 
constructions and programming languages. We hope it will serve as a useful tool to further under- 
stand the nature and structure of quantum computation. The remainder of the paper is organized 
as follows. Section |2 presents the traditional model of quantum computing and its implementation 
in Haskell, focusing on the possibility of structuring the effects using monads. Section [3] discusses 
the limitations of the traditional model as a complete model of quantum computation which should 
include measurement. Section 0] introduces a more general model of quantum based on density 
matrices and superoperators. Our main result is discussed in Section [S] where we show that gen- 
eral quantum computations including measurement can be structured using the generalization of 
monads called arrows. Section gives two complete examples implementing a Toffoli circuit and 
the teleportation experiment: both examples use the arrow notation to express the structure of the 
computation elegantly. Section discusses the limitations of our model and its connection to the 
functional quantum programming language QML 2, . Section[S]concludes. Appendix lAl explains the 
basics of the Haskell notation used in the paper, and the next two appendices present the proofs 
that are omitted from the main body of the paper. 

2 The Traditional Model of Quantum Computing 

We present the traditional model of quantum computing in this section. 
2.1 Vectors 

A finite set a can be represented in Haskell as an instance of the class Basis below. Given such 
a set a representing observable (classical) values, a pure quantum value is a vector a — > C which 
associates each basis element with a complex probability amplitude. The basis elements must be 
distinguishable from each other which explains the constraint Eq a on the type of elements below: 



class Eq a => Basis a where basis : : [a] 
type PA = Complex Double 
type Vec a = a -» PA 



The type constructor Vec is technically not a monad: it corresponds to a Kleisli structure [3] . Yet 
as noted by Mu and Bird (2001), the probabilities introduced by vector spaces constitute a compu- 
tational effect which can be structured using a slight generalization of monads in Haskell 9_. From 
a programming perspective, a monad is represented using a type constructor for computations m 
and two functions: return :: a — > m a and ~^*= :: m a — > (a — > m b) — > m b. The operation ~^*= 
(pronounced "bind" ) specifies how to sequence computations and return specifies how to terminate 
computations: 



return : 


Basis a => a -> Vec a 




return a 


b = if a=b then 1 else 




(»=) :: 


Basis a => Vec a -» (a -» Vec b) 


-» Vec b 


va »= f 


= A b -» sum [ (va a) * (f a b) 


1 a € basis] 



Because of the additional constraint that our computations must be over specified bases whose 
elements must be comparable, the types of our operations are more restricted than strictly desired 
for a monad. However return and 3= satisfy the three monad laws. 

Proposition 2.1 Vector spaces satisfy the required equations for monads. 

Proof. See Appendix [BJ □ 



Vector spaces have additional properties abstracted in the Haskell class MonadPlus. Instances 
of this class support two additional methods: mzero and mplus which provide a "zero" computation 
and an operation to "add" computations: 



mzero : : 


Vec a 






mzero = 


const 






mplus : : 


Vec a -» 


Vec a -» 


Vec a 


mplus v_ 


1 v_2 a = 


v_l a + 


v_2 a 


mminus : 


: Vec a -» Vec a -» 


Vec a 


mminus v 


_1 v_2 a 


= v_l a 


- v_2 a 



For convenience, it is also possible to define various kinds of products over vectors: the scalar 
product $*, the tensor product (*), and the dot product (•): 



($*) : : PA Vec a -> Vec a 
pa $* v = Aa -» pa * v a 

(<*>) : : Vec a -» Vec b -» Vec (a,b) 
vl <*> v2 = A (a,b) vl a * v2 b 

(<.>) : : Basis a => Vec a -» Vec & -> PA 

vl <.> v2 = sum (map (Aa -» conjugate (vl a) * (v2 a)) basis) 



Examples of vectors over the set of booleans may be defined as follows: 

instance Basis Bool where basis = \_False,Bool\ 
qFalse,qTrue,qFT,qFmT :: Vec Bool 
qFalse = return False 
qTrue = return True 

qFT = (1 / sqrt 2) $* (qFalse 'mplus' qTrue) 
qFmT = (1 / sqrt 2) $* (qFalse 'mminus' qTrue) 

The first two are unit vectors corresponding to basis elements; the last two represent state which are 
in equal superpositions of False and True. In the Dirac notation, these vectors would be respectively 
written as \ False), \ True), -^(\False)+ \True)), and -^(\ False) — \ True)). 

Vectors over several values can be easily described using the tensor product on vectors or the 
Cartesian product on the underlying bases: 

instance (Basis a, Basis b) => Basis(&, b) where 
basis = [(a, b) I a € basis, b € basis ] 



pl,p2,p3,epr :: Vec (Bool, Bool) 



pi = qFT <*> qFalse 
p2 = qFalse <*> qFT 
p3 = qFT <*> qFT 

epr (False, False) = 1 / sqrt 2 
epr (True, True) = 1 / sqrt 2 

In contrast to the first three vectors, the last vector describes an entangled quantum state which 
cannot be separated into the product of independent quantum states. The name of the vector "epr" 
refers to the initials of Einstein, Podolsky, and Rosen who used such a vector in a thought experiment 
to demonstrate some strange consequences of quantum mechanics [S]. 

2.2 Linear Operators 

Given two base sets A and B a linear operator / S A —o B is a function mapping vectors over A to 
vectors over B. We represent such operators as functions mapping values to vectors which is similar 
to representation used by Karczmarczuk (2003): 



type Lin a b = a -> Vec b 

fun21in : : (Basis a, Basis b) => (a -» b) -» Lin a b 
fun21in f a = return (f a) 



The function fun2lin converts a regular function to a linear operator. For example, the quantum 
version of the boolean negation is: 

qnot :: Lin Bool Bool 
qnot = fun21in -i 

Linear operations can also be defined directly, for example: 

phase :: Lin Bool Bool 

phase False = return False 

phase True = (0 :+ 1) $* (return True) 

hadamard : : Lin Bool Bool 
hadamard False = qFT 
hadamard True = qFmT 

The definition of a linear operation specifies its action on one individual element of the basis. 
To apply a linear operation / to a vector v, we use the bind operation to calculate v S= /. For 
example (qFT 3= hadamard) applies the operation hadamard to the vector qFT which one can 
calculate produces the vector qFalse as a result. 

It is possible to write higher-order functions which consume linear operators and produce new 
linear operators. An important example of such functions produces the so-called controlled opera- 
tions: 

controlled :: Basis a => Lin a a -> Lin (Bool, a) (Bool, a.) 

controlled f (bl,b2) = (return bl) <*> (if bl then f b2 else return b2) 

The linear operator / is transformed to a new linear operator controlled by a quantum boolean value. 
The modified operator returns a pair whose first component is the input control value. The second 



input is passed to / only if the control value is true, and is otherwise left unchanged. For example, 
(qFT (*) qFalse) (controlled qnot) applies the familiar controlled-not gate to a vector over two 
values: the control value is a superposition of False and True and the data value is False. As one 
may calculate the result of this application is the epr vector. 

Linear operations can be combined and transformed in several ways which we list below. The 
function )*( produces the linear operator corresponding to the outer product of two vectors. The 
functions linplus and lintens are the functions corresponding to the sum and tensor product on 
vectors. Finally the function o composes two linear operators. 



adjoint : : Lin a b -» Lin b a 
adjoint f b a = conjugate (f a b) 

(>*<) : : Basis a => Vec a -» Vec a -» Lin a a 
(vl >*< v2) al a2 = vl al * conjugate (v2 a2) 

linplus : : (Basis a, Basis b) => Lin a b -» Lin a b -» Lin a b 
linplus f g a = f a 'mplus' g a 

lintens :: (Basis a, Basis b, Basis c, Basis d) => 
Lin a b -» Lin c d -» Lin (a,c) (b,d) 
lintens f g (a,c) = f a <*> g c 

o : : (Basis a, Basis b, Basis c) => Lin a b -» Lin b c -» Lin a c 
o f g a = (f a »= g) 



2.3 Example: A Toffoli Circuit 
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The circuit diagram uses the de-facto standard notation for specifying quantum computations. Each 
line carries one quantum bit (qubit); we refer to the three qubits in the circuit as top, middle, and 
bottom. The values flow from left to right in steps corresponding to the alignment of the boxes which 
represent quantum gates. The gates labeled H, V, VT, and Not represent the quantum operations 
hadamard, phase, adjoint phase, and qnot respectively. Gates connected via a bullet to another wire 
are controlled operations. 

In general all three qubits in the circuit may be entangled and hence the state vector representing 
them cannot be separated into individual state vectors. This means that, despite the appearance to 
the contrary, it is not possible to operate on any of the lines individually. Instead the circuit defines 
a linear operation on the entire state: 

toffoli :: Lin (Bool , Bool , Bool) (Bool , Bool , Bool) 
toffoli (top, middle .bottom) = 
let cnot = controlled qnot 

cphase = controlled phase 
caphase = controlled (adjoint phase) 
in hadamard bottom >>= A bl -» 

cphase (middle, bl) »= A (ml,b2) -» 



cnot (top, ml) »= A (tl,m2) -> 
caphase (m2,b2) »= A (m3,b3) -» 
cnot (tl,m3) »= A (t2,m4) -t 
cphase (t2,b3) »= A (t3,b4) -• 
hadamard b4 »= A b5 -» 
return (t3,m4,b5) 



3 Measurement 

The use of monads to structure the probability effects reveals an elegant underlying structure for 
quantum computations. This structure can be studied in the context of category theory and exploited 
in the design of a calculus for quantum computation (141 1151 1131 15] . 

Unfortunately in the traditional model of quantum computing we have used so far, is difficult 
or impossible to deal formally with another class of quantum effects, including measurements, deco- 
herence, or noise. We first give one example where such effects are critical, and then discuss various 
approaches in the literature on how to deal with such effects. 

3.1 Teleportation 

The idea of teleportation is to disintegrate an object in one place making a perfect replica of it 
somewhere else. Indeed quantum teleportation 0] enables the transmission, using a classical com- 
munication channel, of an unknown quantum state via a previously shared epr pair. 

In the following diagram, Alice and Bob initially have access to one of the qubits of an entangled 
epr pair, and Alice aims to teleport an unknown qubit q to Bob: 
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The calculation proceeds as follows. First Alice interacts with the unknown qubit q and her half 
of the epr state. Then Alice performs a measurement collapsing her quantum state and getting two 
classical bits m\ and m<i that she transmits to Bob using a classical channel of communication. 

Upon receiving the two classical bits of information, Bob interacts with his half of the epr state 
with gates controlled by the classical bits. The circuit in the figure can be shown to re-create the 
quantum state q which existed at Alice's site before the experiment. 

Our main interest in this circuit is that it is naturally expressed using a sequence of operations on 
quantum values which include a non-unitary measurement in the middle. Using the model developed 
in the previous section, it is not possible to describe this algorithm as stated. In the next section, 
we briefly several possible ways to deal with this problem. 



3.2 Dealing with Measurement 



The literature includes several approaches to the problem of measurement. We characterize such 
approaches in three broad categories: deferring measurements, using classical control with pointers 
and side-effects, and using density matrices and superoperators. We discuss the first two approaches 
in the remainder of this section, and expand on the latter approach in the next section. 

3.2.1 Deferring measurements: 

The first approach (used for example by Mu and Bird (2001), Van Tonder (2003; 2004) and Karcz- 
marczuk (2003) relies on the principle of deferred measurement 8 . This principle can be used to 
transform computations to always defer measurements to the end. Using this idea one can focus 
entirely on managing the probability effects and perform the measurements outside the formalism. 
The drawback of this approach is clear: programs that interleave quantum operations with mea- 
surements cannot be expressed naturally. For example, transforming the teleportation circuit above 
to defer the measurements until after Bob's computation completely changes the character of the 
experiment, because no classical information is transmitted form Alice to Bob. 

3.2.2 Classical Control and Side-effects: 

In general, this category of models is based on the so-called QRAM (quantum random access ma- 
chine) model of Knill (1996), which is summarized by the slogan "quantum data, classical con- 
trol" 12; . In this context, a quantum computer can be seen as a classical computer with a quantum 
device attached to it. The classical control sends instructions for the quantum machine to execute 
unitary operations and measurements. A measurement collapses the quantum (probabilistic) com- 
putation and forces it to produce a classical (deterministic) result. In fact, the situation is even more 
complicated: measuring part of a quantum state collapses not only the measured part but any other 
part of the global state with which it is entangled. The most common approach to computationally 
realize this hybrid architecture is via manipulating what arc effectively pointers to a global shared 
quantum state as the following examples show: 

• In the flowchart notation for the language introduced by Selinger (2004), the state is repre- 
sented by a collection of variables that can each be assigned once. An operation can only be 
applied to an initial group of the variables (and is implicitly composed with the identity on 
the remaining variables). If the variables are not in the desired order, they must be permuted 
first. Thus the first few steps of the toffoli circuit are: 



^input ql,q2,q3 : qubit 





ql, ql, q3 : qubit 


permute <|) 1 






q3, q2, ql: qubit 


q3,q2,ql *_ Hxld 






q3, q2, ql: qubit 


permute <j>2 






ql, q3, ql : qubit 


q2, q3, ql *= cV x Id 





• In the procedural language QCL a quantum register is a realized using pointers to the 
complete state. Operations on a register map to operations on the state as follows. If we have 
an m-qubit register r which points to an n-qubit state, then an operation U on the register is 
realized using: 

U{r) = 14 (U x I(n - to)) IT r 

The operation U is composed with the identity on the remaining number of qubits of the state. 
The operator II r is an arbitrary reordering operator and nj is its inverse. After re-ordering, 
the lifted U composed with the identity is applied, and the result is permuted back to the 
original order. 

• Jan Skibihski (2001) produced an early Haskell simulator of a quantum computer. The simu- 
lator maintains quantum registers and allows operations to act on specific qubits using what is 
essentially pointers. To apply an operation to the third, fifth, and seventh qubits on a quantum 
register, some low-level calculations depending on the indices and size of the register are used 
to produce a lifted operation composed with several identity operations that acts on the entire 
register. 

• Valiron et. al. (2004) develop a functional quantum programming language based on the 
original work of Selinger (2004). The representation of quantum data in their calculus uses an 
external n-qubit state Q. Programs may contain free variables which are essentially pointers 
to the quantum state. 

• In our previous work jl 1| we introduced virtual values to hide the management of pointers to 
the global state. Using virtual values the code for the toffoli example is essentially identical to 
the one presented earlier, except for the need to manually generate the adaptors which mediate 
between the virtual value and the global state. 

The use of pointers and sharing to model the side-effect of measurement is in some sense ade- 
quate. However by doing so, we completely lose the monadic structure and the direct connections 
to categorical semantics. 

4 Density Matrices and Superoperators 

Fortunately the usual model of quantum computing can be generalized to solve the problem of 
modeling measurements in a better way. In the generalized model, the state of the computation 
is represented using a density matrix and the operations are represented using superoperators pQ. 
Using these notions, the projections necessary to express measurements become expressible within 
the model. We review this model in this section. 

4.1 Density Matrices 

Intuitively, density matrices can be understood as a statistical perspective of the state vector. In the 
density matrix formalism, a quantum state that used to be modeled by a vector v is now modeled 
by its outer product. 



type Dens a = Vec (a, a) 

pureD : : Basis a => Vec a -» Dens a 
pureD v = lin2vec (v >*< v) 

lin2vec : : (a -» Vec b) -» Vec (a,b) 
lin2vec = uncurry 
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The function pureD embeds a state vector in its density matrix representation. For convenience, 
we uncurry the arguments to the density matrix so that it looks more like a "matrix." For example, 
the density matrices corresponding to the vectors qFalse, qTrue, and qFT can be visually represented 
as follows: 

1/2 1/2 
1/2 1/2 

The appeal of density matrices is that they can represent states other than the pure ones above. 
In particular if we perform a measurement on the state represented by qFT, we should get False 
with probability 1/2 or True with probability 1/2. This information which cannot be expressed 
using vectors, can be represented by the following density matrix: 

1/2 A / \ _ / 1/2 
J + \ 1/2 ) ~ \ 1/2 

Such a density matrix represents a mixed state which corresponds to the sum (and then normal- 
ization) of the density matrices for the two results of the observation. If we further calculate with 
the result of measuring qFT by for example, applying the hadamard operation, we get one of the 
two vectors qFT or qFmT, each with probability 1/2. Because all operations on vectors are linear, 
we can express this step as follows: 

1/2 \ =H ( 1/2 \ ( ^ _ ( 1/2 



1/2 J ^ J 1 \ 1/2 J \ 1/2 

As the calculation shows, the application of hadamard has no effect on the density matrix, and indeed 
there is no observable difference between the two configurations before and after the application of 
hadamard. Indeed, the density matrix representation loses the information in the state vectors that 
is not observable ^2] an d hence is a better representation from a semantic perspective. 



4.2 Super-operators 

Operations mapping density matrices to density matrices are called superoperators: 



type Super a b = (a 


a) 


-» Dens b 






lin2super : : (Basis 


a, 


Basis b) => Lin 


a b -> Super 


a b 


lin2super f (al,a2) 




(f al) <*> (dual 


(adjoint f) 


a2) 


where dual f a 


b 


= f b a 







The function linlsuper constructs a superoperator from a linear operator on vectors. To understand 
the basic idea, consider the density matrix resulting from the application of / to \v). This corresponds 
to the outer product of the vector / \v) with itself, which applies / to \v) and the adjoint of / to 
the "dual vector." 



4.3 Tracing and Measurement 

In contrast to the situation with the traditional model of quantum computing, it is possible to 
define a superoperator which "forgets" , projects, or traces out part of a quantum state as well as a 
superoperator which measures part of a quantum state: 

trL : : (Basis a, Basis b) => Super (a,b) b 

trL ((al,bl) , (a2,b2)) = if al s a2 then return (bl,b2) else mzero 



meas : : Basis a => Super a (a, a) 

meas (al,a2) = if al = a2 then return ( (al , al) , (al , al) ) else mzero 



For example, the sequence: 
pureD qFT >>= meas >>= trL 

first performs a measurement on the pure density matrix representing the vector qFT. This mea- 
surement produces a vector with two components: the first is the resulting collapsed quantum state 
and the second is the classical observed value. The last operation forgets about the collapsed quan- 
tum state and returns the result of the classical measurement. As explained earlier the resulting 
density matrix is: 

f 1/2 \ 
\ 1/2 J 

5 Superoperators as Arrows 

By moving to density matrices and superoperators, it becomes possible to express both the original 
computations as well as measurements in the same formalism. One might hope that the original 
monadic structure of quantum computations is preserved, but it appears that this is not the case. 
The best we can do is to prove that the new model of computation fits within a generalization of 
monads called arrows. 



5.1 Arrows 

The application of a superoperator to a density matrix can still be achieved with the monadic bind 
operation, instantiated to the following type: 

>>= : : Dens a -» ((a, a) -» Dens b) -» Dens b 



This type does not however correspond to the required type as computations now consume 
multiple input values. This observation is reminiscent of Hughes's motivation for generalizing monads 
to arrows 7 . Indeed, in addition to defining a notion of procedure which may perform computational 
effects, arrows may have a static component independent of the input, or may accept more than one 
input. 

In Haskell, the arrow interface is defined using the following class declaration: 

class Arrow a where 

arr :: (b-»c)-»abc 

(»>) :: abc->acd->abd 

first :: a b c -» a (b,d) (c,d) 

In other words, to be an arrow, a type a must support the three operations arr, ^>, and first with 
the given types. The operations must satisfy the following equations: 



arr id 3g> / = / 

/ ^> arr id = f 

(f^g)^h = f^>(g^> h) 

arr (<?./) = arr / ^> arr g 

first (arr /) = arr (f x id) 

first (f^g) = first f^> first g 

first f ^> arr (id x g) = arr (id x g) ^> first f 

first f ^> arr fst = arrfst ^> / 

first (first f) ^> arr assoc — arr assoc 3S> first f 



where the functions x and assoc are defined as follows: 



(/ x g) (a,b) = (/ a, g b) 
assoc ((a, b), c) = (a, (6, c)) 

Graphically the functions associated with the arrow type are the following: 
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The function arr allows us to introduce "pure" arrows which are simple functions from their 
inputs to their outputs. The function ;g> is similar to it composes two computations. The 
function first is the critical one for our purposes: it allows us to apply an arrow to a component of 
the global quantum state. The equations above ensure that these operations are always well-defined 
even with arbitrary permutations and change of associativity. 



5.2 Superoperators are Arrows (with Eq constraint) 

Just as the probability effect associated with vectors is not strictly a monad because of the Basis 
constraint, the type Super is not strictly an arrow as the following types include the additional 
constraint requiring the elements to be comparable: 

arr : : (Basis b, Basis c) => (b -» c) -» Super b c 
arr f = fun21in (A (bl,b2) -» (f bl, f b2)) 

(>>>) :: (Basis b, Basis c, Basis d) => 

Super b c -» Super c d -» Super b d 

(»>) = o 



first :: (Basis b, Basis c, Basis d) => Super b c -» Super (b,d) (c,d) 
first f ((bl.dl) , (b2,d2)) = permute ((f (bl,b2)) <*> (return (dl,d2))) 
where permute v ( (bl ,b2) , (dl ,d2) ) = v ( (bl ,dl) , (b2 ,d2) ) 

The function arr constructs a superoperator from a pure function by applying the function to both 
the vector and its dual. The composition of arrows is simply the composition of linear operators. 
The function first applies the superoperator / to the first component (and its dual) and leaves the 
second component unchanged. The definition calculates each part separately and then permutes the 
results to match the required type. 

Proposition 5.1 Superoperators satisfy the required equations for arrows. 

Proof. See Appendix [O D 



The proposition implies that we can use the arrow combinators to structure our computations. 
For instance, the first few steps of the Toffoli circuit of Section would now look like: 



toffoli : : Super (Bool , Bool , Bool) (Bool , Bool , Bool) 
toffoli = 

let hadS = lin2super hadamard 

cphaseS = lin2super (controlled phase) 
cnotS = lin2super (controlled qnot) 
in arr (A (aO, bO, cO) -» (cO, (aO, bO))) »> 

(first hadS >» arr (A (cl, (aO, bO)) -» ((bO, cl) , aO))) »> 
(first cphaseS »> arr (A ((bl, c2) , aO) -» ((aO, bl) , c2))) »> 
(first cnotS »> arr (A ((al, b2) , c2) -» ((b2, c2) , al))) »> ... 

Clearly this notation is awkward as it forces us to explicitly manipulate the entire state and to 
manually permute the values. However, all the tedious code can be generated automatically as we 
explain next. 

5.3 A Better Notation for Arrows 

Following the Haskell's monadic do-notation, Paterson (2001) presented an extension to Haskell with 
an improved syntax for writing computations using arrows. We concentrate only on the explanation 
of new forms which we use in our examples. Here is a simple example to illustrate the notation: 

el : : Super (Bool, a.) (Bool, a.) 
el = proc (a,b) -» do 

r <— lin2 super hadamard -< a 

returnA -< (r,b) 

The do-notation simply sequences the actions in its body. The function returnA is the equivalent 
for arrows of the monadic function return. The two additional keywords are: 

• the arrow abstraction proc which constructs an arrow instead of a regular function. 

• the arrow application -< which feeds the value of an expression into an arrow. 

Paterson (2001) shows that the above notation is general enough to express arrow computations 
and implemented a preprocessor which translates the new syntax to regular Haskell. In the case of 
el above, the translation to Haskell produces the following code: 

e2 : : Super (Boo I, a.) (Boo I, a) 
e2 = first (lin2super hadamard) 

As the example shows, the output of the preprocessor is quite optimized. 

5.4 Superoperators are (probably) not monads 

Arrows are more general than monads. In particular, they include notions of computation that 
consume multiple inputs as well as computations with static components, independent of the input. 
Due to this general aspect of arrows, there are some subclasses of them which turns out to be 
equivalent to monads. More precisely, arrow types which support the following app function are just 
as expressive as monads. 

class Arrow => ArrowApply a where 
app : : a (a b c, b) c 

In other words, for superoperators to be monads, we would have to define a superoperator of type: 
Super (Super be, b) c 

which in our case would require Super b c to be an instance of Basis. Unfortunately there is no 
straightforward way to view the space of superoperators as a finite set of observables. 



6 Examples Revisited: Toffoli and Teleportation 



Using arrows and the notation introduced by Patterson, we can express both of our examples ele- 
gantly. 

6.1 Toffoli 

The code mirrors the structure of the circuit and the structure of the monadic computation expressed 
earlier: 

toffoli :: Super {Bool , Bool , Bool) (Bool , Bool , Bool) 
toffoli = let hadS = lin2super hadamard 

cnotS = lin2super (controlled qnot) 
cphaseS = lin2super (controlled phase) 
caphaseS = lin2super (controlled (adjoint phase)) 
in proc (aO,bO,cO) -» do 
cl <— hadS -< cO 
(bl,c2) <— cphaseS -< (bO.cl) 
(al,b2) <— cnotS -< (aO,bl) 
(b3,c3) <— caphaseS -< (b2,c2) 
(a2,b4) <— cnotS -< (al,b3) 
(a3,c4) <— cphaseS -< (a2,c3) 
c5 <— hadS -<; c4 
returnA -< (a3,b4,c5) 

6.2 Teleportation 

We use the machinery we have developed to faithfully express the circuit presented in Section |3. II 
We break the algorithm in two individual procedures, alice and bob. Besides the use of the arrows 
notation to express the action of superoperators on specific qubits, we incorporate the measurement 
in Alice's procedure, and trace out the irrelevant qubits from the answer returned by Bob. 

alice :: Super (Bool, Bool) (Bool, Bool) 
alice = proc (eprL.q) -» do 

(ql,el) <— (lin2super (controlled qnot)) -< (q.eprL) 

q2 <— (lin2super hadamard) -< ql 

((q3,e2) , (ml,m2)) <— meas -< (q2,el) 

(ml',m2') <- trL ( (q3 , e2) , (ml ,m2) ) 

returnA -< (ml' ,m2') 

bob :: Super (Bool , Bool , Bool) Bool 
bob = proc (eprR,ml,m2) -• do 

(m2',el) <— (lin2super (controlled qnot)) -< (m2,eprR) 

(ml',e2) <— (lin2super (controlled z)) -< (ml, el) 

q' <- trL -< ((ml' ,m2') ,e2) 

returnA -< q' 

teleport :: Super (Bool , Bool , Bool) Bool 
teleport = proc (eprL,eprR,q) -• do 

(ml,m2) «— alice -< (eprL.q) 

q' <— bob -< (eprR,ml,m2) 

returnA -< q' 



7 Linear Typing: QML 



The category of superoperators is considered to be an adequate model of non-reversible quantum 
computation J2| ■ O ur construction presented so far seems to suggest that this category corresponds 
to a functional language with arrows, and so that we can accurately express quantum computation 
in such a framework. But as we explain below, this is not quite the whole story. 

First consider the well-known "non-cloning" property of quantum states |3| . The arrow notation 
allows us to reuse variables more than once, and we are free to define the following operator: 

copy :: Super Bool (Bool, Bool) 
copy = arr (A x -» (x,x)) 

But can this superoperator be used to clone a qubit? The answer, as explained in Section 1.3.5 of 
the classic book on quantum computing 8 , is no. The superoperator copy can be used to copy 
classical information encoded in quantum data, but when applied to an arbitrary quantum state, 
for example like qFT, the superoperator does not make two copies of the state qFT but rather it 
produces the epr state which is the correct and desired behavior. Thus, in this aspect the semantics 
of arrows is coherent with quantum computation, i.e., the use of variables more than once models 
sharing, not cloning. 

In contrast, in our model there is nothing to prevent the definition of: 

weaken :: Super (Bool, Bool) Bool 
weaken = arr (A (x,y) -» y) 

This operator is however not physically realizable. Applying weaken to epr gives qFT. Physically 
forgetting about x corresponds to a measurement: if we measure the left qubit of epr we should get 
qFalse or qTrue or the mixed state of both measurements, but never qFT . 

Therefore, our use of Haskell as a vehicle for expressing the ideas finally hits a major obstacle: 
arrow computations must be required to use every value that is introduced. Instead of attempting to 
continue working within Haskell, a better approach might be to now consider a functional quantum 
language like QML whose type system is designed to explicitly control weakening and decoherence, 
and to express the separation of values and arrow computations in that framework. 

In more detail, QML is a functional quantum programming language which addresses this 
problem by using a type system based on strict linear logic: contraction is permitted and modelled 
by copy while weakening has to be explicit and is translated by a partial trace. QML also features to 
case operators: a classical case operator which measures a qbit and returns the appropriate branch 
and a quantum case operator which avoids measurement but requires that the branches return results 
in orthogonal subspaces. 

QML programs can be compiled to quantum circuits, using the category of finite quantum 
computation FQC — Grattage's QML compiler [H] is based on this semantics. An irreversible 
computation can be modelled by a reversible circuit, allowing additional heap qubits, which are 
initialized to a predefined values at the beginning of the computation and disposing, i.e. measuring, 
qbits at the end of the computation. To any FQC morphism we can assign a superoperator and 
indeed every superoperator can be represented this way. 

Alternatively, we can interpret QML programs directly as superoperators, giving rise to a con- 
structive denotational semantics exploiting the library of arrow combinators developed here. We 
hope to exploit this semantics to further analyze QML and to develop high level reasoning principles 
for QML programs. 

8 Conclusion 

We have argued that a realistic model for quantum computations should accommodate both unitary 
operations and measurements, and we have shown that such general quantum computations can 



be modeled using arrows. This is an extension of the previous-known observation that one can 
model pure quantum probabilities using monads. Establishing such connections between quantum 
computations and monads and arrows enables elegant embeddings in current classical languages, 
and exposes connections to well-understood concepts from the semantics of (classical) programming 
languages. We have demonstrated the use of arrows to model elegantly two examples in Haskell, 
including the teleportation experiment which interleaves measurements with unitary operations. 
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A A Haskell Primer 

We use Haskell as a precise mathematical (and executable) notation. 

It is useful to think of a Haskell type as representing a mathematical set. Haskell includes several 
built-in types that we use: the type Boolean whose only two elements are False and True; the type 
Complex Double whose elements are complex numbers written a :+ b where both a and b are elements 
of the type Double which approximates the real numbers. Given two types a and 6, the type (a, b) 
is the type of ordered pairs whose elements are of the respective types; the type a — > b is the type 
of functions mapping elements of a to elements of b; and the type [a] is the type of sequences (lists) 
whose elements are of type a. For convenience, we often use the keyword type to introduce a new 
type abbreviation. For example: 

type PA = Complex Double 

introduces the new type PA as an abbreviation of the more verbose Complex Double. A family of 
types that supports related operations can be grouped in a Haskell class. Individual types can then 
be made an instance of the class, and arbitrary code can require that a certain type be a member 
of a given class. 

The syntax of Haskell expressions is usually self-explanatory except perhaps for the following 
points. A function can be written in at least two ways. Both the following definitions define a 
function which squares its argument: 

sq n = n * n 

sq' = A n -» n * n 

A function / can be applied to every element of a list using map or using list comprehensions. If xs 
is the list [1, 2, 3, 4], then both the following: 

map sq xs 

[ sq x I x <— xs ] 

evaluate to [1,4, 9, 16]. 

Usually, a function / is applied to an argument a, by writing / a. If the function expects two 
arguments, it can either be applied to both at once / (a, b) or one at a time fab depending on its 
type. When convenient the function symbol can be placed between the arguments using back quotes 
a '/' b. 



B Proof of Monad Laws for Vectors. 

Proof of Proposition ^. II The definitions of return and satisfy the three monad laws: 

• First monad law: (return x) ^= / = / x 

(return x) ^$*= f = X b . sum [ return x a * / a b \ a <— basis] 

= X b . sum [ if x = a then 1 else * / a b \ a <— basis] 
= Xb.f xb 
= fx 

• Second monad law: m ^= return = m 

m return — X b . sum [ma* return a b \ a <— basis] 

= A b . sum [ma * if a = b then 1 else | a <— basis] 
= X b . m b 

= m 

• Third monad law: (m 3= /) ~S*= g = m 3*= (A x . f x 3= g) 

(m /) ^= g = (Xb . sum [ma * f ab[ a «— basis]) ^= g 

= X c . sum[(sum [ma * f a b \ a <— basis]) * g b c \ 
b <— &asis] 

= A c . sum [m a * f a b * g b c \ a <— basis, b <— 6asis] 

m (X x . f x ^3= 9) = Ac. sum[m a * (/ a ^3= 3) c | a <— 6asis] 

= A c .swm[m a * (sum[f 5 * g b c \ b <— 6asis]) 
a <— 6asis] 

= A c .sum [ma * f a b * g 6 c | a <— basis, b <— 6asis] 

C Proof of Arrow Laws for Superoperators 

Proof of Proposition 

• First arrow equation: arr id ^> f = f ■ 

arr id ^> / = fun2lin (X (al, a2).(id al, id a2)) 'o'/ (by arr and ^>) 
= fun2lin id 'o' / (6j/ simplification) 

— return ' 0' / (by fun2lin) 
= X a . return a ^3= / (6j/ '0') 

= X a . f a (by monad law 1.) 

- / 

• Second arrow equation: / ^> arr id = f . 

f ^> arr id = f V fun2lin (X (61, 62) . (id 61, id 62)) (6y arr and 

= / ' o' fun2lin id (by simplification) 

= f l o l return (by fun2lin) 

= X a . f a ^3= return (by o) 

= X a . f a (by monad law 2.) 
= f 



• Third arrow equation: (/ ^> g) ^> ft = f ^> (g ft). 

(/ »> g) »> h = (/ V ff ) V ft (fey »>) 

= Xb . (Xa . f a »= g) b »= ft (by o) 
= A6.(/6 >^ .9) >^ ft (by (3) 

f »> (5 » h ) = / 'o' ( 5 '0' ft) (6» >H 

= Xa.f a »= (Xb . g b »= h) (by 0) 
= X a . (f a 3= .9) ^= ft (6j/ monad law 3.) 

• Fourth arrow equation: arr (g . f) = arr f arr g. 

arr (g . f) = fun2lin (X (61, 62) .((# . /) 61, ( 5 . /) 62)) (6t/ arr) 

= return .(A (61, 62) . ((3 . /) 61, (g . f) 62)) (by fun2lin) 
= X (61, 62) . return ((g . /) 61, (5 . /) 62) (simplification) 

arr J ^> arr .g = /«n2fe (A (61, 62) . (/ 61, / 62)) V fun2lin(X (61, 62) . (.9 61, 5 62)) 
(6y ^3= and arr) 

= return . (X (61, 62) . (/ 61, / 62)) V return .(A (61, 62) . (.9 61, g 62)) 
(by fun2lin) 

= X (61, 62) . retern (/ 61, / 62) »= A (61, 62) . return (g 61, <? 62)) 
(6y 0) 

= A (61, 62) . (A (61, 62) . return (g 61, g 62)) (/ 61, / 62) 

(by monad law 1.) 
= A (61, 62) . return ((g . f) 61, (.9 . /) 62) (by (i) 

• Fifth arrow equation: first (arr /) = arr (f x id). 

first (arrf) = first (fun2lin (X (61, 62) . (/ 61, / 62))) (by arr) 

= first (return . (X (61, 62) . (/ 61, / 62))) (by fun2lin) 

= first (X (61, 62) . return (/ 61, / 62)) (by simplification) 

= X ((61, dl), (62, d2)) . X ((x, y), (w, z)) .return (f 61, / 62) (x, w) * 

return (dl , d2) (y , z) (by first) 
= X ((61, dl), (62, d2)) . X ((x, y), (w, z)) . 

if ((/ 61,/ 62), (dl, d2)) == ((a, w), (2/, 2)) then 1 else (by return) 

arr (f x id) = = fun2lin (X ((61, dl), (62, d2)). ((/ 61, dl), (/ 62, d2))) (by arr) 

= return . (X ((61, dl), (62, d2)). ((/ 61, dl), (/ 62, d2))) (by fun2lin) 
= X ((61, dl), (62, d2)) . return ((/ 61, dl), (/ 62, d2)) 
= X((bl,dl),(b2,d2)).X((x,y),(w,z)). 

if ((/ 61, dl), (/ 62, d2)) == ((x, y), (w, zj) then 1 else (by return) 

• Sixth arrow equation: first (f ^> g) = first f first g. In the following proofs assume: 
adl ((61,dl),(62, d2)) = (61,62) and ad2 ((61, dl), (62, d2)) = (dl, d2) . 

first (j 'o' g) = first (X a . f a »= g) 

(by l o l ) 

= A 6 . A ((x, y), (w, zj) . (f (adl b) ^= g) (x, w) * return (ad2 6) (y, z) 
(by first) 

= A 6 . A ((x, y), (w, z)) .(Ac. sum [(/ (adl 6)) a * g a c \ a <— 6asis]) (x, w) 

* return (ad2 6) (y, z) (by >=) 
= A 6 . A ((x, y), (w, zj) . sum[(f (adl bj) a * g a (x,w) \ a <— 6asis] * 

return (ad2 6) (y, z) (by (3) 



first f'o' first g = A a . first f a ^= A b . first g b 
(by 'o') 

= A a . A ((x, y), (w, z)) . f (adl a) (x, w) * return (ad2 a) (y, z) 3= 
A b . A ((a;, y), (w, z)) . g (adl b) (x, w) * return (ad2 6) (y, z) 
(by first) 

= A a . A ((x, y), (w, z)) . sum [f (adl a) (m, o) * return (ad2 a) (n,p)* 

(A ((a;, y), (w, z)) . g (m, o) (x, w) * return (n,p) (y, zj)((x, y), (w, z)) 

((m,n),(o,p)) <— basis] (by »=) 
= A a . A ((a;, y), (w, z)) . sum [ / (adl a) (m,o) * return (ad2 a) (n,p)* 

g(m,o)(x,w) * return (n, p) (y, z) \{{m, n), (o, p)) <— basis] 
= A a . A ((x, y), (w, z)) . sum [/ (adl a) al * g al (x, w) * 

return (ad2 a) a2 * return a2 (y, z) \ al <— basis , a2 <— basis] 

(by simplification) 

= A a . A ((x, y), (w, z)) . sum [/ (adl a) al * g al (x, w) \ al <— basis ] 
* return (ad2 a) (y, z) (by simplification) 

• Seventh arrow equation: first f ^> arr (id x g) = arr (id x g) ^> first f . 

Ihs = first f 'o' arr (id x g) 

lhs= A ((al, 61), (a2, b2)) . first f ((al, bl), (a2, 62)) 3^ 

fun2lin (A ((a, 6), (c, d)) . ((a, g 6), (c, g d))) (6y 'o' and arr) 
= A ((al, 61), (a2, 62)) . first f ((al, 61), (a2, 62)) ^ 

A ((a, 6), (c, d)) . return ((a, g 6), (c, 5 d)) (by fun2lin) 
= A ((al, 61), (a2, 62)) . A ((x, y), (w, z)) . f (al, a2) (x, w) * return (bl, 62) (y, z) »= 

A ((a, b), (c, d)) . return ((a, g b), (c, g d)) (by first) 
= A ((al, bl), (a2, b2)) .Ac. sum [/ (al, a2) (m, o) * return (bl, 62) (n, p) * 

return ((m, gn),(o, g p)) c\((m,n),(o,p)) <— basis] (by 
= A ((al, bl), (a2, b2)) . A ((#, y), (w, z)) . sum [f (al, a2) (m, 0) * return (61,62) (n,p) * 

return ((m, g n),(o, g pj) ((x,y),(w,z))\ ((m,n),(o,pj) <- basis] 

(by simplification) 
= A ((al, 61), (a2, 62)) . A ((x, y), (w, z)) . sum [f (al, a2) (m, o) * 

[if (61, 62) == (n, p) then 1 else 0] * 

[(if (m, g n), (o, g p)) == ((x, y), (w, z)) then 1 else 0] | ((m, n), (o,p)) <— basis] 
(by return) 

= A ((al, 61), (o2, 62)) .A ((x, y), (w, z)) . if (g bl, g 62) == (y, z) 
then / (al, a2) (x, w) else 

rhs = arr (id x g) 'o' /irsi / 

rhs = A ((al, 61), (a2, 62)) . fun2Un (A ((a, 6), (c, d)) . ((a, g b), (c, g d))) 
((al, bl), (a2, 62)) ^ first f (by V and arr) 
= A ((al, 61), (a2, 62)) . return ((al, g bl), (a2, g 62)) ^ first f 
(by Jun2lin) 

= A ((al, bl), (a2, b2)) . first f ((al, g bl), (a2, g b2)) (by monad law 1.) 
= A ((al, bl), (o2, 62)) .A ((x, y), (w, z)) . f (al, o2) (x, w) * return (g bl, g 62) (y, 2) 
(by /irst) 

= A ((al, bl), (a2, 62)) .A ((x, y), (w, z)) . f (al, o2) (s, w) * 
[if (5 61, g 62) == (y, 2) then 1 elseO] (by return) 



• Eighth arrow equation: first f ^> arr fst = arr fst ^> f. 

Ihs = first f 'o' arr(X(a, b).a) 

Ihs = A((al,61),(a2, 62)) . first f ({al,bl),(a2, 62)) »= arr A (a, 6) . a (6y o) 

= A ((al, 61), (a2, 62)) . first f ((al, 61), (a2, 62)) »= A ((a, 6), (c, d)) . return (a, c) 
(6y arr ) 

= A ((al, 61), (a2, 62)) . A ((x, y), (w, z)) . f (al, a2) (x, w) * 

return (61, 62) (y, z) ^= A ((a, 6), (c, d)) . return (a, c) (by first) 

= A ((al, 61), (a2, 62)) . A (cl, c2) . sum [/ (al, a2) (m, o) * return (61, 62) (n,p)* 
return (m, o) (cl, c2) | ((m,n),(o,p)) <— basis] (by ^=) 

= A ((al, 61), (a2, 62)) . A (cl, c2) . sum [f (al, a2) (m, o) * 
[if (61, 62) == (n,p) then 1 else 0] * 

[if (m, o) == (cl, c2) then 1 else 0] | ((m, n), (o, p)) <— basis] (by return) 
= X ((al, 61), (a2, 62)) . A (cl, c2) . / (al, a2) (cl, c2) (by simplification) 

rhs = arr fst l o l f 

rhs = A ((a, 6), (c, d)) . return (a, c) V / (by arr) 

= A ((al, 61), (a2, 62)) . (A ((a, 6), (c, d)) . return (a, c)) ((al, 61), (a2, 62)) ^ f 
(by o) 

= A ((al, 61), (a2, 62)) . / (al, a2) (by monad law 1.) 
= A ((al, 61), (a2, 62)) . A (cl, c2) . / (al, a2) (cl, c2) 

• Ninth arrow equation: first (first f) ^> arr assoc = arr ass oc ^> first f 

Ihs = A(((al, 61), cl), ((a2, b2),c2)).f Hrst( first /)(((al, 61), cl), ((a2, 62), c2)) 3^ 
arr(A((a,6),c).(a, (6, c))) 



Ihs = A (((ol, 61), cl), ((a2, 62), c2)) . first (X b . \ ((x, y), (w, z)) ./ (adl 6) (x, w) * 
return (ad2 6) (y,z)) (((al, 61), cl), ((a2, 62), c2)) »= 
A (((al, 61), cl), ((a2, 62), c2)) . return ((al, (61, cl)), (a2, (62, c2))) 
(6j/ first) 

= A (((al, 61), cl),((a2, 62), c2)) . A ((ml,nl),pl) ((m2, n2), p2) . 

(A 6 . A ((:r, ?/), (w, 2)) . / (adl 6) (x, w) * return (ad2 6) (y, z)) ((al, 61), (a2, 62)) 
((ml,nl),(m2,n2)) * return (cl, c2) (pl,p2) »= A (((al, 61), cl), ((a2, 62), c2)) . 
return ((al, (61, cl)), (a2, (62, c2))) {by first) 

= A(((al,61),cl),((a2, 62), c2)) . A ((ml, nl),pl) ((m2, n2), p2) . 

/ (al, a2) (ml, m2) * return (61, 62) (nl, n2) * return (cl, c2) (pl,p2) »= 
A (((al, 61), cl), ((a2, 62), c2)) . return ((al, (61, cl)), (a2, (62, c2))) 

(byP) 

= A (((al, 61), cl), ((a2, 62), c2)) . A ((xl, (yl, zl)), (x2, (y2, z2))) . 

sum [f (al, a2) (ml, m2) * return (61 , 62) (nl,n2) * return (cl, c2) (pi, p2) * 
retern ((ml, nl),pl) ((m2, n2), p2) ((xl, (yl, zl)), (a;2, (y2, z2))) 
((ml, nl), pi) ((m2, n2), p2) <— basis] 
(by 3-) 

= A (((al, 61), cl), ((o2, 62), c2)) . A ((xl, {yl, zl)), (x2, (y2, z2))) . 
sum [/ (al, a2) (ml, m2) * [if (61, 62) == (nl, n2) then 1 else 0] * 
[if(cl,c2) == (pl,p2)then 1 else 0] * 

[if ((ml, nl), pi) ((m2,n2), p2) == {{xl, {yl, zl)), {x2, (y2, z2))) then 1 else 0]| 
((ml, nl), pi) ((m2, n2), p2) <— basis] 
(by return) 

= A (((al, 61), cl), ((a2, 62), c2)) . A ((xl, (yl, zl)), (x2, (y2, z2))) . f (al, o2) (xl, x2) * 
retern ((61, cl), (62, c2)) ((yl, zl), (y2, z2)) 



rhs = X(((al, 61), cl), ((a2, 62), c2)) . return((al, (61, cl)), (a2, (62, c2))) 'o' first f 

rhs = A(((al,61),cl),((a2, 62), c2)) . return ((al, (61, cl)), (a2, (62, c2))) »= ^rsi/ 
(6y 0) 

= A (((al, 61), cl), ((a2, 62), c2)) . first f ((al, (61, cl)), (a2, (62, c2))) 

(by monad law 1.) 
= A (((al, 61), cl), ((o2, 62), c2)) . A ((xl, (yl, zl)), (x2, (y2, z2))) . 

f (al, o2) (xl, x2) * retern ((61, cl), (62, c2)) ((yl, zl), (y2, z2)) (by first) 
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